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Thin cylindrical tethers are common lipid bilayer membrane structures, arising in situations rang- 
ing from micromanipulation experiments on artificial vesicles to the dynamic structure of the Golgi 
apparatus. We study the shape and formation of a tether in terms of the classical soap-film problem, 
which is applied to the case of a membrane disk under tension subject to a point force. A tether 
forms from the elastic boundary layer near the point of application of the force, for sufficiently large 
displacement. Analytic results for various aspects of the membrane shape are given. 
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I. INTRODUCTION 

Imagine a soap film connecting two nearby parallel 
rings which are slowly pulled apart. The film evolves 
through a series of catenoids until the ring separation 
reaches a critical value. Then the film breaks. Now imag- 
ine performing the same experiment with microscopic 
rings connected by a lipid bilayer membrane. In this 
paper we will show that the membrane forms catenoidal 
shapes at small ring separations, but instead of break- 
ing, the membrane forms a thin cylindrical tether for 
sufficiently large displacement. 

A situation very much like this thought experiment 
arises in a host of real experiments on artificial vesicles, 
living cells, and organelle membranes such as that of the 
Golgi apparatus. Perhaps the most controlled tether ex- 
periment is that of Evans and Yeung, in which a tether 
forms when a micropipet is withdrawn from a spherical 
vesicle held at a fixed tension and bonded to a stationary 
bead 0. There are many variations on this experimen- 
tal theme ||^, |^, ^, |^, |^ , and such experiments have been 
used to measure a wide variety of membrane mechanical 
properties 0. 

Tethers commonly form in less controlled situations as 
well. Simply pulling on a vesicle or cell with a sufficiently 
large point force leads to a membrane tether (Fig. |^). 
Tethers form when tubulin trapped inside a vesicle poly- 
merizes to form microtubules ||): as the microtubules 
grow, the initially spherical membrane at first distorts 
into an ellipsoidal shape, and then eventually forms a 
surface of revolution with a contour in the shape of the 
Greek letter "0" [|[ |ll|. Improvements in staining 
techniques have recently revealed dynamic tether net- 
works in the Golgi apparatus of living cells fl^ ; similar 
model membrane networks have been studied in vitro and 
used as templates for making more durable networks p3[ . 



The variety of conditions under which tethers occur 
shows that they are robust structures, insensitive to the 
details of applied forces and boundary conditions. De- 
spite much theoretical, computational, and experimental 
attention, a simple intuitive picture for tether formation 
which exploits the insensitivity to geometric constraints 
(such as volume conservation) and emphasizes the role of 
the most essential control parameters has yet to emerge. 
It is therefore natural to study the theory of tethers in the 
simplest possible situation, namely the classic soap film 
geometry described above (Fig. ^) . This model situation 
captures the essential features of the tether shape and for- 
mation without the experimentally important, but ulti- 
mately complicating, effects of volume and lipid conserva- 
tion. Our quasi-analytic approach exploits the smallness 
of the tether radius, and is complementary to important 
numerical work by various groups, most notably Hein- 
rich et al. We begin our analysis in section II with a 
discussion of shells, balloons, and soap films in order to 
contrast their familiar mechanical properties with the pe- 
culiar properties of lipid bilayer membranes. The latter 
are the subject of section III, which defines precisely the 
model problem. Section IV reviews the elements of the 
soap film problem which are relevant for understanding 
the formation of tethers, the subject of section V. There 
we use asymptotic methods to solve the linearized equa- 
tions for small ring displacements and large tensions, and 
find that the membrane shape is that of a catenoid with 
a small elastic boundary layer surrounding the smaller 
ring (the point force). At a critical ring separation, a 
tether forms from the elastic boundary layer. The tether 
shape is studied analytically and numerically, and a con- 
nection is made with the classic film-coating calculations 
of Landau and Levich iQ, as well as Bretherton's re- 
lated calculation of the shape of a long air bubble rising 
in a fluid-filled capillary tube jl^. Section VI. is the con- 
clusion. The appendix reviews the subtleties that arise 
when comparing the variational approach for lipid bilayer 
membrane elasticity to that of moment and force balance. 
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FIG. 1: Equilibrium shapes of a vesicle subject to various point forces. An optical tweezer exerts the point force, which increases 
from left to right. Figure courtesy of D. Fygenson. 



II. SHELLS, BALLOONS, AND SOAP FILMS 



Before exploring the lipid bilayer membrane properties 
that give rise to tethers, we review three examples of elas- 
tic surfaces encountered at the macroscopic scale: shells, 
balloons, and soap films. Much intuition about elastic 
surfaces derives from these canonical examples. Some of 
this intuition can be directly applied to lipid membranes, 
but it is illuminating to point out the crucial differences 
as well. 

Shells are solid surfaces with a small thickness and a 
preferred shape in the absence of external stresses. A 
plate is a shell with a flat preferred shape. It is a matter 
of common experience that the force required to bend a 
thin plate through a certain displacement is much less 
than that required to stretch the plate through the same 
displacement. The plate can undergo large displace- 
ments through bending without subjecting each element 
to large stress; thus, linear elasticity is valid and the non- 
linearities in the equations for shape are solely geometri- 
cal. In the limit of small plate thickness, stretching en- 
ergy exceeds bending energy, since bending is differential 
stretching. Deformations without stretching are there- 
fore of lowest energy fl^ . These deformations, in which 
the distances between nearby points remain fixed, are 
called isometric. For example, the axisymmetric isomet- 
ric deformations of the plane are cylinders and cones [ p7| . 
Note that a shearing motion locally stretches the plate, 
even if the total area does not change. Bending energy 
(and possibly boundary conditions) removes the degen- 
eracy when, as in the case of a flat plate, there is a mul- 
tiplicity of isometric deformations. Bending also plays 
a role near boundaries and point forces. A familiar ex- 
ample of a plate with vanishing thickness is a sheet of 
paper, which bends easily but hardly stretches or shears. 
The reader can easily verify that a sheet of paper sub- 
ject to a point force forms a cone with a single fold, a 
shape with two-fold symmetry. For shells, the isometric 
constraint on deformations is even more severe. Any de- 
formation of a sphere requires some local stretching or 



shearing: no almost-spherical shapes are isometric to the 
sphere. A spherical shell is geometrically rigid. On the 
other hand, there are isometric small deformations of a 
spherical shell with a circular hole The problem 

of determining the geometric rigidity of an axisymmet- 
ric shell with boundaries was solved using the qualitative 
theory of differential equations iQ; for recent progress 
using a geometric approach, which can be generalized to 
non-axisymmctric shells, see ||l9| . 

Next, consider the case of balloons. Like plates with 
vanishing thickness, there is virtually no cost to bend- 
ing a section of a balloon compared to that of stretch- 
ing or shearing. However, deformations are not required 
to be isometric since the cost of stretching is also low. 
Such a material is called a "membrane" in the theory 
of elasticity iQ; to avoid confusion with lipid bilayer 
membranes, this terminology will not be used. Since 
real balloons are easily stretched out of the linear elas- 
tic regime, the nonlinearity of the governing equations 
arises from the constitutive relations as well as the ge- 
ometry of large deformations. Anyone who has inflated 
a cylindrical balloon has seen a phenomenon very much 
like the tethers of Fig. ^ A partially inflated balloon 
has two cylindrical regions, one smaller than the other, 
which are smoothly connected by a junction region. How- 
ever, this phenomenon differs fundamentally from that of 
tether formation in lipid bilayer membranes for several 
reasons. No point force is required to make a balloon 
form a tether; it arises from the nonlinear constitutive 
relation between tension and areal extension, which is 
reflected by a non-monotonic p-V curve, reminiscent of 
the isotherms of the van der Waals equation of state for a 
gas . The material in the larger cylinder is stretched 
more, and thus has a higher tension. The smaller cylin- 
der is under smaller tension; the difference between the 
axial force of each tension is precisely balanced by the 
net pressure on the junction region. Thus, a pressure 
jump across the surface of the balloon is crucial for this 
tether. Below it is shown that such a pressure jump is 
unnecessary for lipid bilayer membrane tethers. 

Our final canonical example is the "ideal" soap film. 
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FIG. 2: Model problem. Two parallel rings with aligned cen- 
ters are connected by a lipid bilayer membrane in contact with 
a lipid reservoir at fixed chemical potential (per unit area) /i. 



which simply minimizes area subject to the boundary 
conditions and volume constraints. Effects such as thick- 
ness variations and draining do not have a clear coun- 
terpart in the case of lipid bilayer membranes, and are 
therefore disregarded. Since soap films are liquid, the 
static in-plane shear rigidity vanishes, and the notion of 
geometric rigidity does not apply. Deformations need 
not be isometric and folds are not required to bend a 
flat film into a hemisphere. (There can be solid-like fold- 
ing behavior in liquid film dynamics whenever bending 
flows are preferable to extensional flows [0.) Unlike the 
case of plates, interfacial tension is a material property of 
the soap film, and is not determined by external forces. 
To imagine poking a soap film with a point force, con- 
sider again the geometry of two rings described in the 
introduction, but now take one ring radius to be much 
smaller than the other. Equilibrium solutions exist only 
when the ring separation is less than a critical separa- 
tion, comparable to the radius of the smaller ring. At a 
slightly larger separation, the film breaks. The equilib- 
rium shapes are shallow catenoids, and look nothing like 
tethers (structures reminiscent of thin tethers form dur- 
ing the rupture of a soap film, but these have a dynam- 
ical origin; see However, these shallow catenoidal 
shapes resemble the junction region between the tether 
and vesicle of Fig. ^ Soap films and catenoids will be 
studied more fully in section IV. 

Lipid bilayer membranes have elements of each of these 
examples. Like macroscopic plates, lipid membranes re- 
sist bending, but they are typically fluid and therefore 
lack geometric rigidity, like soap films. The interplay be- 
tween these solid-like and fluid properties leads to the 
formation of tethers. 



III. LIPID BILAYER MEMBRANES 
A. Bending elasticity 

There is a vast literature on the mechanics of lipid bi- 
layer membranes {e.g. see p^). Here we review only 
the aspects relevant to tethers. Lipid molecules are am- 
phiphilic, composed of two oily hydrophobic chains at- 
tached to a polar hydrophilic head. Lipids self-assemble 
in aqueous solution to shield the chains from the water, 
forming micron-size bilayer surfaces. In the fluid phase, 
a typical self-diffusion constant for a lipid molecule in a 
membrane is of order 10~®cm^/sec The membrane 
is clearly a two-dimensional fluid over the time scale of 
typical tether experiments: a molecule will diffuse from 
one side of a one-micron radius spherical vesicle to the 
other in about a second. Hence, there is no in-plane 
static shear modulus. However, membranes have a non- 
zero bending modulus since bending the membrane com- 
presses and extends the slightly elastic heads and tails of 
the lipid molecules. The fluid nature of the membrane 
constrains the form of the elastic energy to depend only 
on the shape of the membrane. To lowest order in cur- 
vatures, the bending energy is given by the expression of 
Canham and Helfrich 



26 



dS 



(1) 



where H is the mean curvature of the membrane surface, 
and K is the Gaussian curvature (k and k are elastic 
constants). In terms of the principal radii of curvature 
R\ and i?2 at a point on the surface, IH = l/i?i + 
l/i?2 and K = l/(i?ii?2)- Explicit formulas for these 
curvatures will be given below. For simplicity, suppose 
that there is no difference between the two sides of the 
bilayer; hence, there is no spontaneous curvature. The 
Gaussian curvature term is typically dropped in studies 
of vesicles since it amounts to the sum of a deformation- 
independent term and a boundary term by the Gauss- 
Bonnet theorem: 



dSK = 47r(2 ~2g)+ (b ds n. 



(2) 



where the genus g is the number of handles and Kg is the 
geodesic curvature of the boundary of the surface [ p7[ . 
Note that a more detailed treatment of the formation of 
tethers in vesicles would require additional terms of the 
generalized bilayer couple model |2^, |3^, In 
the spirit of explaining tether formation in the simplest 
possible context, these are disregarded. 

Although typical values for the elastic bending mod- 
ulus K are lO-lSfc^r, thermal fluctuations easily excite 
long-wavelength bending modes since the elastic energy 
vanishes as the fourth power of the wavenumber for fluc- 
tuations about a flat sheet. These fluctuations lead to 
an entropic area elasticity similar to that of semiflexible 
polymers, as is most directly illustrated by the experi- 
ments of Evans and Rawicz [E2j (see also p3|). In these 
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experiments, the tension in a vesicle is measured as a 
function of apparent area by suctioning a small amount 
of the vesicle membrane into a pipet. At low suctions (or 
low tensions) , the resistance to stretching is the entropic 
penalty of reducing the number of fluctuating modes. 
At high tensions, most of the thermal ripples have been 
smoothed out, and the resistance to stretching is mainly 
due to the membrane's intrinsic area elasticity. Thermal 
fluctuations are therefore negligible for the high tension 
regime considered here. 



(recall i? = 1), the problem is to minimize the energy 
(measured in units of /ii?^) 

£ = j dS+^j dS{2Hf (4) 

for given ring separation /iq subject to the boundary con- 
ditions and e ^ 1. Equation (^ casts the tether problem 
into the same form as the classic variational problem for 
a minimal surface (e = 0). 



B. Model problem and nondimensionalization 



C. Euler-Lagrange equations 



As in the introduction, consider a lipid membrane 
spanning two rings that are initially concentric and ly- 
ing in the plane z = (see Fig. |^). The outer ring of 
radius R remains in this plane, but the height of the 
inner ring will be varied. Measuring all lengths in units 
of the radius _R, henceforth R = 1. The inner ring then 
has radius rp <C 1. Later, we will take ro — > to model 
the application of a point force. For comparison with the 
case of soap films and also for the numerical approach, it 
is convenient to keep tq nonzero for now. Assume that 
there is a reservoir of lipid at a fixed chemical potential 
(per unit area) /i. Further suppose this "surface ten- 
sion" is very large compared to the bending elasticity, 
IJlR^ ^ K. Thermal fluctuations are therefore negligible. 
Since the membrane has edges, the Gaussian modulus R 
affects the shape through the boundary conditions (recall 
the geodesic curvature term of Eq. (y)). For simplicity, 
we disregard this effect and set k = 0. 

Our model is closest in spirit to the tether experiments 
of Evans and Yeung [Q mentioned earlier. The pipet suc- 
tion sets the value of the tension /i, and the small amount 
of lipid projecting inside the pipet serves as a reservoir. 
In the experiments, the vesicle is under pressure and is 
therefore curved, whereas in our model there is no pres- 
sure jump across the initially flat membrane. The role 
of the pressure in vesicle tethers is small since the tether 
curvature is much larger than the vesicle curvature for 
high tension. In fact, it is shown in section VB2 that the 
pressure leads to a subleading correction to the tether 
radius. 

The complete specification of the model problem re- 
quires defining the boundary conditions at the rings. The 
tangent plane of the surface at the point force is perpen- 
dicular to the force if the line of action is along the axis 
of symmetry. Since the role of the small ring is to mimic 
a point force, the membrane is clamped so that the tan- 
gent plane at each point on the boundary (with the small 
ring) is in the plane of the small ring. The ring at r = 1 
is somewhat artificial; therefore we choose the simplest 
possible boundary condition for r = 1, which turns out 
to be zero moment, H = The outer ring acts as a 
hinge. Defining the dimcnsionless parameter 



(3) 



The derivation of the Euler-Lagrange equations from 
the energy of Eq. (^ is somewhat lengthy but straight- 
forward |34[|. The result is 



2e{V^H + 2H'^ - 2HK) - 2i? + Ap = 0, 



(5) 



where is the covariant Laplacian on the surface. Ap, 
measured in units of /i/i? (in other words, /i) is zero for 
our soap film geometry, but it is included in Eq. (||) for 
later discussion of the effect of pressure on tether shape. 
Note that e = and Ap = yields the minimal-surface 
equation, H — 0. Since the membrane shape is a surface 
of revolution, natural coordinates for the surface are ip, 
the azimuthal angle, and s, arclength along a meridian. 
Arclength is measured from the inner ring, which has 
coordinate s = 0. The position of a point on the surface 
is therefore X(s, tp) = r(s)f -I- z{s)z, where r and z are 
cylindrical coordinates. Note that + = 1, since s is 
arclength. 

With these choices, the metric, or first fundamental 
form, is 

gijdCd^^ =ds^ + r^d(p'^, (6) 
where = s and ^■^ = ip. The second fundamental form 



Kijd^d^^ = {zsTss - rsZss) ds^ - rzsdip'^ 



(7) 



where Zg — dz/ds, etc. We follow the usual conventions 
for raising and lowering indices using the inverse g^^ of 
the metric tensor. Thus, g^^gkj — <5!, = g^'^K^j, and 
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1 d d 

r ds ds ' 



(8) 
(9) 
(10) 



where g is the determinant of the metric tensor gij , and 

+ z1 = \ was used to simplify (0) . 

The Euler-Lagrange equation ^) for an axisymmet- 
ric shape is a nonlinear ordinary differential equation, 
easily solved with numerical methods. However, more 
insight is gained by exploiting the smallness of e. Since 
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e multiplies the term in with the highest number of 
derivatives, there will be a boundary layer. The layer 
occurs where the boundary conditions are incompatible 
with Eq. (|^) with e = 0. It is now clear why the zero 
moment H = boundary condition is the most natural 
condition at the larger ring: this boundary condition is 
compatible with the minimal-surface equation and does 
not lead to a boundary layer at the larger ring. The 
clamped boundary condition at the smaller ring, how- 
ever, is incompatible with the minimal-surface equation, 
since the ring must exert a moment on the surface to 
keep it clamped. In this boundary layer, bending and 
tension balance and the shape is determined by the full 
Euler-Lagrange equation, where the smallness of e is off- 
set either by large s-derivatives or large curvatures. Bal- 
ancing eW^H with H in (1q) reveals that the boundary 
layer thickness scales as ^/e. In the outer region, beyond 
the elastic boundary layer, bending is unimportant and 
the shape is governed by the minimal-surface equation, 
H = 0. The only nonplanar axisymmetric minimal sur- 
face is a catenoid, the surface of revolution generated by 
a catenary p7| . Thus, the membrane forms a catenoid in 
the outer region. In the next section, reviews basic facts 
about catenoids. 



IV. CATENOID LORE 

We noted in the previous section that setting e = 
and Ap = in Eq. results in the minimal-surface 
equation H — 0. To leading order, the outer solution 
to Eq. (||) is given by exactly the same condition. The 
solution to the minimal-surface equation is conveniently 
found by applying Noether's theorem directly to the en- 
ergy functional (H) with e — 0. To this end, rewrite dS* 
in terms of r{z): 




1.0 b 



FIG. 3: Ring separation ho vs. minimum neck radius b (force 
by 2ir) for the values vq = 0.2,0.4,0.6,0.8, and 1.0, proceed- 
ing from the innermost curve to the outermost curve. For a 
given ro, as the rings are separated, the sequence of shapes 
corresponds to the trajectory marked with arrows. 



The boundary condition 
± cosh" ^(1/6), and thus Eq. 



r{0) = 1 determines c 
(p^) expands to 



cosh ^■^^ -F \/l — 6^ sinh 



(14) 



The upper sign is chosen as it corresponds to catenaries 
with a minimum neck radius at a positive value of z (c > 
0). The force F, or equivalently the minimum neck radius 
b, is determined by the boundary condition at the other 
ring: r{ho) = ro, or 



cosh 



\/ 1 — sinh 



(15) 



Solving (|T^) for ho gives the separation as a function of 
force: 



S — dipdz r-y/l -I- ' 



(11) 



ho^b log 



^0 ± 



fe2 



1- 



(16) 



The conserved quantity associated with the invariance of 
the integrand of ( [ll] ) with respect to translations in z is 
the axial force F necessary to hold the rings apart at a 
given separation: 



F 
2^ 



Vi+^ 



(12) 



The radius r attains its minimum value b — F/ {2tt) when 
Tz = 0. Integrating ( |l2|) yields the catenoid 



6 cosh 



z - 



c 



(13) 



where c is the ^-coordinate of the minimum radius. Note 
that b is the minimum possible radius; i.e. it is possibly 
not attained. The minimum radius is attained only if 
c < ho- 



Note that the two branches form a closed curve in the 
b-ho plane for tq < 1 (Fig. ^). Since each curve has a 
maximum (marked with a dot for the curves tq — 0.6 
and ro = 1.0), there is a critical ro-dependent separa- 
tion beyond which no catenoidal solution exists. The 
soap film spanning the two rings breaks just beyond this 
critical separation. For fixed ro and a given separation 
ho below the maximum, there are two catenoidal solu- 
tions. For example, Fig. ^illustrates the two equilibrium 
catenoids with ro = 1 and ho = 0.6. For a given ho, one 
can show that the catenoid with the larger b has less area. 
Thus, the solution with the smaller neck {e.g. the upper 
catenoid in Fig. ^) is unobservable in real soap films. In 
the presence of bending stiffness and a fully-developed 
tether, it is shown below that the axial force is 27r-\/2e. 
To see which catenoid matches onto a fully developed 
tether, consider the extreme case rg = 1. Since the ax- 
ial force vanishes as e — > 0, the matching catenoid must 
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with 



FIG. 4: The two catenoidal solutions with equal-size rings 
(ro = 1) and ring spacing ho — 0.6; the lower catenoid has 
less area. 



have b near zero and lie on the left branch of the vq — 1 
curve of Fig. By continuity, as rg decreases, the corre- 
sponding matching catenoid is always on the left branch. 
Therefore, the catenoid joining a tether is the one with 
a narrow neck; bending stiffness selects this otherwise 
unobservable shape. 

Finally, the arrows in Fig. ^ display the trajectory 
of shapes as the rings are separated. The force F in- 
creases from zero to a maximum value, and then de- 
creases slightly before the film breaks. This nonmono- 
tonic behavior occurs in the case of tethers as well, as is 
shown below. 



V. FROM CATENOIDS TO TETHERS 

As discussed earlier, for small but nonzero e, the outer 
region of the membrane will have a catenoidal shape, and 
there will be an elastic boundary layer near the small 
ring. This elastic boundary layer allows the attainment 
of the point force limit, rp — s- 0, in contrast to the case 
of the soap film. We begin with the analysis of the case 
of small axial separation, ho <^ 1. 



A. Small displacements 

When ho <^ 1, it is most convenient to work in the 
Monge parameterization, in which the surface is repre- 
sented by its height z(r) above the plane z = Q. To lead- 
ing order in ho, s « r, the mean curvature H w (1/2) V^z, 
and the Euler-Lagrange equation reduces to 



0, 



(17) 



r dr 



dr 



(18) 



The boundary conditions at the large ring 



1 



and the condition of zero moment, V z = 0. 



are 
At 



the inner ring r = vq, the displacement z = ho and the 
slope dh/dr — 0. 

At order e°, the outer solution satisfies V^Zoutor(?') — 0, 



I.e. 



ZoutcA'^) =bi+ 62 log r. 



(19) 



The boundary condition on displacement at r = 1 fixes 
61 = 0. The zero- moment boundary condition adds no 
constraint on the solution of (pj|); 62 must be determined 
by matching to the inner solution. 

Note that the outer solution diverges at the inner ring 
as ro — > 0; the inner solution must correct for this di- 
vergence. To find the inner solution, expand the region 
near r — with the rescaling p = rj^fk. Then the inner 
solution satisfies 



1 d 
p dp \^ dp 



pAp\ dp) 



0, (20) 



Ci +C2logp + C3/o(p) + C4i^o(p), (21) 



where /o(p) and Ko{p) are modified Bessel functions. C3 
must vanish, since /o(p) diverges as p — > 00 and cannot be 
matched to the outer solution. The boundary conditions 
at r = rg (i.e., p = rp/v^) add two more constraints to 
yield 



= 1- 



C4 



Koip)-Ko 

THurA— 



(22) 



To match the inner and outer solutions, note that 
Ko{p) decays exponentially at large p. Therefore, the 



constant terms of z;. 



must vanish, and the coefficient 



of the logarithmic term of Zj, 
ing order. 



must match 62- To lead- 



■^-outer 

ho 



^ifi(^)logr 



^ifi(^)logro + Xo(^)^ 



and 



ho ^^i(^)logro + i^o(^) 



(23) 



(24) 



The Bessel function Ko{r/y/e) cancels the logarithmic di- 
vergence of the log(r) term of the inner solution. To con- 
struct a uniformly- valid approximation Zcompositc('') for 
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FIG. 5: Comparison of membrane profiles 2compositc(r) for 
ro = 0.1 (upper curve) and the point force limit ro 
(lower curve). 



both the inner and outer regions, add the two solutions 
and subtract their common part This procedure 

yields the very compact result 



^'Composite 



(25) 



It is now clear why a small amount of bending elasticity 
allows a membrane to support a point force: the elastic 
boundary cuts off the divergence of the logarithm of the 
outer solution. To find the solutions in the limit of a 
point force, recall i^i(ro/^/e) = Ve/j'o + C(^o/Ve) and 
^o(ro/V^) = -7 + log(2Vi/ro) + O((ro/V^)2), where 
7 = 0.5772... is the Euler constant. Thus 



^composite (^) — ^0 



logr + Xo(r/^) 
-7 + log(2Vi) 



(26) 



Note that the boundary condition 2compositc(l) = is 
satisfied up to terms of order e^^'* exp(— l/y^)/ loge for 
small e. Figure || shows that although the details of the 
membrane shape depends somewhat sensitively on the 
value of ro, the physical limit — > is well behaved. In 
this limit, the elastic boundary layer becomes a small disc 
of approximate radius \ft around the point force. Thus, 
the outer solution is roughly the catenoid that connects 
a ring of radius unity with a ring of radius -^/e. As in the 
previous section, the maximum ring separation for a soap 
film in this situation is approximately equal to the radius 
of the smaller ring. Thus, as /iq increases, the amplitude 
of the catenoid increases until is of order y^. Since 
the amplitude of the catenoid cannot increase beyond this 
value, the boundary layer deforms into a thin cylinder to 
accommodate further increases in h^. The formation of 
the tether is a smooth process; there is no bifurcation. 



B. Tether: analytical approach 

1. Tether radius, tether stability, and axial force 

Tether formation is an intrinsically nonlinear phe- 
nomenon, and to give a complete account of the tether 
shape we resort to numerical methods. However, many 



features of the tether yield to an analytic approach. The 
tether radius is the most prominent such feature. Our nu- 
merical calculations will verify that the tether has a cylin- 
drical shape between from the end cap and the catenoidal 
junction. Thus, for our soap-film geometry with Ap = 0, 
the radius follows from (||) with constant mean curvature 
and vanishing Gaussian curvature: 



H ^0. 



(27) 



Since H = — l/(2a) for a cylinder of radius a, the exact 
tether radius a — y/e/2 ||l| . This square-root dependence 
of tether radius on inverse tension has been verified ex- 
perimentally by Evans and Yeung [Q . 

To study the stability of a tether, write r{z) = a + u{z) 
and expand the elastic energy (^) to O(u^) (to express 
the metric (^ and mean curvature (^) as functions of z 
use dz/ds = ^1 + r'{z)^), which yields 



2n 



= d 



e 

2a 



1 - 



202) 



ea 



."2 



,2 \ 1 



(28) 



A total derivative term has been dropped in (|2^). Min- 
imizing the M-independent terms over a yields the equi- 
librium tether radius. The terms linear in u vanish as 
expected when a takes the equilibrium value ■\/e/2. Note 
that the terms quadratic in u' vanish in equilibrium as 
well, since the terms u and au''^ /2 always enter in the 



2U /2. This combination arises from the 
in the original energy. Since the remain- 



combination u- 
f actor r\/l + r' 
ing terms of Eq. ( |28| ) in u'^ and u"^ are positive definite, 
the tether is stable. Therefore, the equilibrium cylinder 
solution does not undergo a pearling instability |3^, ^7\ . 
These considerations suggest that the pearling behavior 
induced by a rapid pull of a vesicle tether arises be- 
cause hydrodynamic resistance prevents the radius from 
instantly assuming the value appropriate to the new value 
of tension. This mechanism differs in detail from that of 
the laser-tweezer-induced instability of membrane tubes 
with fixed volume ^t\ . 

The tether radius determines the axial force on the 
catenoid. The argument hinges on the z-independence of 
the axial force. For an undistorted cylinder with u = 
and equilibrium radius a = \/ e/2, the total energy per 
unit length EjL — 27r\/2e, as easily follows from (|2^). 
Therefore, the axial force F/(27r) — ^/2e saturates to 
a constant value independent of tether length once the 
tether has formed. Since the axial force is independent 
of z, the junction connecting the tether to the ring is a 
catenoid with b = V2e. Note that the minimum attain- 
able radius of the limiting catenoid is twice the radius of 
the tether. Thus, the catenoid cannot smoothly join onto 
the cylindrical tether and there must be a transition re- 
gion (see Fig. Before analyzing this transition region 
in more detail, we consider the role of pressure. 
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z 

FIG. 6: Overlay of a fully developed tether [A) with radius 
a — sj e/2 and a catenoid (5) with minimum radius h = 
\plt. In this example, e = 0.005 and the transition region lies 
roughly between z — 0.2 and z = 0.4. 

2. The effects of pressure are subleading 

Evans and Yeung argued that the pressure jump Ap, 
present in the case of a closed vesicle under tension, 
plays little direct role in determining the tether radius . 
Since a sphere of radius Rq has constant mean curvature 
I/Rq and Gaussian curvature the Euler-Lagrange 

equation ^ in the spherical region of the vesicle reduces 
to the Young-Laplace law, 2H = Ap (even in the presence 
of bending resistance). Measuring lengths in units of Rq 
(for this paragraph only), the Euler-Lagrange equation 
in the region of the tether becomes 



since the pressure jump is everywhere uniform. For small 
e, the three solutions are 

a = ±^/7/2 + e/2 + 0{e^/^) (30) 
a = l/2-e + 0{e^). (31) 

The solution a = + 0{e) is unphysical. The so- 

lution a ~ ■\/e/2 + 0(e) corresponds to the tether in the 
case Ap = 0. Thus, to leading order, the tether radius is 
unchanged and the effect of pressure appears at order e. 
Finally, the solution with a radius near 1/2 corresponds 
to a balance of pressure and tension, and is not relevant 
for tethers. 

Returning to our model problem with Ap = 0, an ap- 
parent paradox arises. Since an axial force is required 
to pull the tether out of the membrane disc, there must 
be a tension in the membrane. This tension is isotropic, 
since the membrane is fluid. But consider the cylindrical 




FIG. 7: Axial force due to the greater circumference of the 
outer leaf. 

portion of the tether between two fixed values of z. If 
this cylinder is cut in half along the long axis (Fig. |^), 
then each apparently experiences a resultant force due 
to the tension. What force balances this tension force if 

The paradox is most readily resolved by comparing 
the Euler-Lagrange equations of the variational approach 
with the equations which follow from force and moment 
balance on a membrane element, given the lipid bilayer 
membrane constitutive relation |l[]. The appendix reca- 
pitulates the comparison between the two approaches to 
the equilibrium shape equations. There it is shown that 
the coefficient /i of the area term in the variational energy 
is the tension only for minimal surfaces: /i = S -I- e-ff^, 
where E is the tension. Since the pressure jump Ap = 0, 
the cylindrical region of the membrane is in a state of 
pure bending, E = 0. But due to the differential stretch- 
ing inherent in bending, the outer sheet is stretched and 
the inner sheet is compressed (Fig. Since the sheets 
are fluid, the tension or compression in each sheet is 
isotropic. The compressive and tensile forces along the 
lines of longitude of the cylinder cancel (consistent with 
Ap — 0). However, since the outer sheet is longer than 
the inner sheet along a line of latitude, there is a net 
axial tension. The axial force at zero pressure jump is 
a manifestation of the liquid properties of lipid bilayer 
membranes. 



3. The tether always necks twice 

We have seen that the shape of a membrane subject to 
a point force and under high tension is best described as a 
boundary-layer problem, with tension dominating in the 
outer region and bending dominating in the inner region. 
This situation is reminiscent of the coating problem stud- 
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ied by Landau and Levich [Q , and of Bretherton's bub- 
ble problem , in which the shapes of different regions 
of an interface are determined by different balances. The 
analogy goes further: Bretherton showed that the trail- 
ing edge of a large air bubble, rising in a capillary tube 
filled with viscous liquid, has a slight ripple ||l5|. It will 
now be shown that there are slight ripples in the shape 
of a lipid membrane at both ends of the cylindrical tether 
region. These ripples have been noticed in the numerical 
work of ref. 1^. 

Since bending and tension are equally important in the 
junction region, we must rescale the variables to balance 
these two effects. If Si is the arclength corresponding to 
a point in the transition region, it is enou gh to assume 
that the radius r{si) is small and close to e/2, without 
any further specification of si. It is therefore natural to 
rescale the radius as in section V., r = Pv^- The further 
rescalings a — {s — si)/y/e and C = (z — z{si))/y/e lead 
to a balance of the bending and tension terms: 




where 



2H-^ - 2HK - = 0, 



(32) FIG. 8: Membrane shape for various ring separations; e 
0.005 and ro = 0.005. 



H = 

k = 



Paa 
P 



l_d_ _d_ 
p dcT^dfT ' 



(33) 
(34) 

(35) 
(36) 



and + Q = 1. Therefore, the transition region is 
governed by the full nonlinear Euler-Lagrange equation, 
and there are no further simplifications arising from the 
smallness of e. However, one can use perturbation the- 
ory to study the shape of the transition region near the 
tether. Let p = 1/^/2 + rj, with 77 ^ 1. To leading order 
in r], Eq. (B3) becomes 



(37) 



Note that Eq. also follows immediately from Eq. ( |2^ ) 
with appropriate rescalings. There are four independent 
solutions to Eq. (|37|), each of the form rja = C'a exp(ip(T), 
where p — ±{l±i) and a = 1, ...,4. The shape near ei- 
ther end of the tether region is an exponentially-damped 
sinusoid with wavelength 2ti/^ and decay length ^/e. 



C. Tether: numerical solution 

The last section shows that a description of the mem- 
brane shape in the junction region requires the solution 
of a nonlinear differential equation with no small param- 
eters, despite the smallness of e. Rather than solve this 
equation numerically and match the solution onto the 
tether and catenoidal regions, we simply solve for the 



complete shape numerically. Standard relaxation tech- 
niques 1^ are used to solve for the shape as a func- 
tion of ring displacement /iq, with a small ring of radius 
To = 0.001 mimicking the point force. Figure ^ displays 
the membrane shape for various /iq. For small ho, the 
shape is well approximated by the linearized catenoid 
with an elastic boundary layer at small radius (see section 
v.). As increases, the amplitude of the catenoid in- 
creases until the limiting catenoid with b — -y/e is reached. 
For larger separations, a tether forms. The axial force as 
a function of displacement is shown in Fig. |9[ Note that 
the force increases to a maximum and then decreases 
slightly before saturating to V2e- This behavior is re- 
flected in Fig. H, where the limiting catenoid lies inside 
the catenoids with slightly lower values of /iq , since these 
catenoids have slightly larger values of the minimum neck 
radius b = F/{27r). Figure ^ shows the ripple in the 
junction region. The radial scale has been magnified for 
clarity. 



VI. CONCLUSION 

We have seen that tethers in our model problem are a 
type of boundary-layer phenomenon. In the cylindrical 
tether region, bending dominates, whereas tension dom- 
inates at larger radii. These insights carry over to the 
more complicated problem of tether formation in closed 
lipid bilayer membrane vesicles, where the quantitative 
details of the force vs. extension will be different, since 
tension depends on extension. An important generaliza- 
tion of the problem considered here would be to study 
membranes with varying degrees of in-plane order, rang- 
ing from liquid-crystalline to solid-like, since the liquid 
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FIG. 9: Force vs. displacement; e = 0.005, ro = 0.005. 



FIG. 11: Forces and moments acting on an element of an 
axisymmetric membrane. 
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FIG. 10: Ripple profile, e = 0.0002. 



nature of fluid membranes is crucial for tether formation. 
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APPENDIX A: PLATE THEORY VS. 
VARIATIONAL PRINCIPLE 

This appendix reviews the force and moment balance 
relations for axisymmetric shells ||l|, E3, and the consti- 
tutive relations for fluid membranes KT The approach is 
equivalent to the variational approach taken in the text, 
and elucidates the apparent paradox discussed in section 
VB2. Figure |ll| shows the forces and moments acting 
on a small element of fluid membrane. Only the forces 
and moments that enter the shape equations are shown. 

is the force per unit length parallel to the meridian 
acting on an element edge along the azimuthal direction. 

is the force per unit length in the azimuthal direction 
acting on an element edge along a meridian. The shear- 
ing force Qm acts along the surface normal on an element 
edge along the azimuthal direction. The external stresses 
Pn and pt are forces per area acting on the element in 
the normal and meridional directions, respectively. The 
curvature along the meridian is Cm — d6'/ds, and the 



curvature in the azimuthal direction is 



suiO/r. 



The balance of forces and moments is just the same as 
in shells. Normal stress balance requires 



1 d 

r ds 



(rQrn)- 



Tangential stresses balance when 



Id Tu,dr 

-Pt = --rv'^rn) 3- + CjnQr. 

r as r as 



(Al) 



(A2) 



where dr/ds = cos 6. Moment balance about the ip-axis 
relates the shearing force Qm to the moments per unit 
length M„i and M^: 



l_d 

r ds 



(rM„ 



My dr 
r ds 



(A3) 
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It is useful to have an expression for the total axial 
force acting on a circle of latitude. Consider the resultant 
of the normal and axial stresses along the axial direction: 



r{pn cos 9— pt sint 



— ( rTm sin 9 -rQ„^ cos 6* ) , (A4) 



where we have used Eqs. (A1,A2). In our model prob- 
lem, the external stresses vanish, Pn — Pt = 0. Thus, 
rTm sin9 — rQ„i cos 9 is a constant, the total axial force by 
27r (compare with the axial force on a soap film, Eq. (E2|)): 



F 
2^ 



(A5) 



Returning to the derivation of the shape equations, 
consider now the tangential forces per unit length Tm 
and Tip. If a; denotes the coordinate across the thickness 
of the membrane, then 



(A6) 
(A7) 



where = / dx, f^ = jT^ dx, Mm = J xt^ dx, and 
= J XT^ dx. Part of the "tension" in the membrane 
comes from the bending moments. 

The constitutive relation for the fluid membrane com- 
pletes the specification of the shape equations. Since the 



fluid nature implies isotropy, = t^- Define the com- 



mon value of tension as t 



T,„. Likewise, 



Mm = M^ = KC, 



(A8) 



where c = Cm + c^. Thus, the shearing force is known 
once the curvature of the membrane is known: 



Qm — ^ 

ds 



Tangential force balance, Eq. (A2), becomes 

d / 1 2^ 



(A9) 



(AlO) 



In the absence of flow, the external tangential stresses 
vanish, and t = f -\- kc^/2 is constant. However, f and c 
need not separately be constant. 

Inserting the constitutive relations into Eq. ( |AlD , the 
normal stress balance becomes 



- ^ -, ^2 1 d / dc 

Pn=TC- -KC{Cm-C^) - '^--TX'^^ 

^ r ds \ ds 



(All) 



But since c = 2H and K = CmC^, Eq. ( All ) reduces to 
the Euler-Lagrange equation (^ with ji = t. 
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